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Abstract 

We compute the 0(aba s ) two-loop corrections to the neutral Higgs boson masses in the 
Minimal Supersymmetric Standard Model, using the effective potential approach. Such cor- 
rections can be important in the region of parameter space corresponding to tan ,8 > 1 and 
sizeable \i. In spite of the formal analogy with the 0(a t a s ) corrections, there are important 
differences, since the dominant effects are controlled by the sbottom-Higgs scalar couplings. 
We propose a convenient renormalization scheme that avoids unphysically large threshold 
effects associated with the bottom mass, and absorbs the bulk of the 0(ai,a s + abat) cor- 
rections into the one-loop expression. We give general explicit formulae for the 0{ai ] a s ) 
corrections to the neutral Higgs boson mass matrix. We also discuss the importance of the 
O(a^) corrections and derive a formula for their contribution to in a simple limiting case. 



1 Introduction 



The existence of a light CP-even neutral Higgs boson is a crucial prediction of the Minimal 
Supersymmetric extension of the Standard Model, or MSSM, and has been one of the most 
active areas of theoretical investigations in the last decade. At the tree level, the masses of the 
neutral CP-even Higgs bosons of the MSSM can be computed in terms of three input parameters: 
the mass m A of the neutral CP-odd particle, the mass m z of the weak neutral gauge boson, and 
the ratio of Higgs vacuum expectation values tan/3 = vijv\ (for a review and references, see 
e.g. (Hi). For tan/3 <C mt/mb, the dominant one-loop corrections are the 0(ctt) ones, where 
at = hf/ (47r) and ht is the superpotential top coupling. Such coupling controls both the top- 
Higgs Yukawa couplings and a number of cubic and quartic stop-Higgs scalar couplings, and leads 
to significant contributions from both top and stop loops 0|- The 0(ab) one-loop corrections 
associated with the superpotential bottom coupling hb, where a& = /;|/(47r), can be numerically 
non-negligible only for tan /3 > 1 and sizeable values of the \x parameter. At the classical level 
hb/ht = (m,b/mt) tan/3, thus we need tan/3 3> 1 to have a& ~ at in spite of mj, <C mj. Moreover, 
and in contrast with the top-stop case, numerically relevant contributions can only come from 
sbottom loops: those coming from bottom loops are always suppressed by the small value of 
the bottom mass. A sizeable value of /i is then required to have sizeable sbottom-Higgs scalar 
interactions in the large tan/3 limit. 

We are now at the stage where the most important genuine two-loop corrections are being 
evaluated: general results have been obtained both for the 0(ata s ) || ||, || and for the 0[pq) 
|3], H corrections. In this paper we move one step further, computing the 0(aba s ) corrections 
and discussing the 0{a\ v ) and 0(atab) ones. For convenience, we evaluate two-loop effects 
directly in the physically relevant limit of large tan /?: 



where is the Fermi constant. As a result, we obtain extremely compact analytical formulae. 
Keeping v\ ^ would only generate more complicated expressions, without adding any relevant 
information. 

The plan of the paper is the following. We first give the analytical result at 0(aba s ) and in 
the DR scheme. We then identify a convenient renormalization scheme that avoids unphysically 
large threshold effects and absorbs the largest 0(aba s + o^at) corrections into the one-loop 
expressions. In particular, we discuss how to use the experimental information on the bottom 
mass, which receives large threshold corrections ||, to extract the value of the renormalized 
coupling hb- We finally present numerical results for some representative parameter choices, and 
conclude with an explicit formula for the 0{aba s + a%) corrections to in a simple limiting 
case. 

2 General formulae and DR results 

The momentum-independent part of the one-loop 0(ab) and two-loop 0{aba s ) corrections to 
the neutral CP-even Higgs boson mass matrix can be obtained by taking the second derivatives 
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of the effective potential [] at its minimum, or by performing appropriate substitutions and limits 
in the <D(ata s ) results of |J|. In the limit of Eq. (|T|), we find: 



(AMI 



cff 

11 

eff 



(AXi) 12 



i 2 2 
n b s 29 b 



1 l2 2 

2 n b s 29 b 

h" 2 2 



pie + F 2i 



(2) 
(3) 
(4) 



Before explaining the meaning of the different symbols, we recall that an important simplification 
occurs if we look at the lightest Higgs eigenvalue, rrih, in the limit m A ^> m z , since in that limit 
Am| coincides with (AvM^Dlf ■ Our conventions are such that, at the classical level, the top 
and bottom quark masses are given by m t = h t V2/V2 and m& = h\,v\/\/2, where the Yukawa 
couplings (ht,hb) and the VEVs (vi,V2) are all taken to be real and positive. In addition, we 
assume [i and A b to be real, but we do not make any assumption on their sign, whereas we 
choose the gluino mass rrig to be real and positive. At the classical level, the sbottom mixing 
angle S20 b = sin 2d~ h is given by 
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> m~ are the two eigenvalues of the 
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where the arrow denotes the large tan/3 limit, and m? 

sbottom mass matrix. Finally, the superscripts in the functions (F, G) indicate the order of the 
loop contribution. At one loop, and in the large tan (3 limit, the only relevant function is 
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where N c = 3 is a color factor. Notice that F u is negative definite. 

We first present our results for F 2 ^ and G 2i in the DR scheme. In other words, we assume 
that the 0{a b ) one-loop contribution is written entirely in terms of DR parameters (masses and 
couplings), evaluated at a certain renormalization scale Q. In units of g 2 Cp N c / '(16 7r 2 ) 2 , where 
C F = 4/3, we find: 
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x The effective potential for vanishing CP-odd fields was computed in ||. To make contact with the physical 
m A , the effective potential should be computed as a function of both CP-even and CP-odd fields, as in M . 
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where the hats on F and G denote DR-quantities, and, here and in the following, (1 

the interchange m~ - "'~ 
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2) denotes 



m- . Notice that, in our limit, the bottom quark only contributes through 



bottom-sbottom-gluino diagrams. The above way of presenting the results is convenient for 
analysing models that predict, via the MSSM renormalization group equations, the low-energy 
DR values of the MSSM input parameters in terms of a more restricted set of parameters, assigned 
as boundary conditions at some scale much larger than the weak scale. One of the parameters, 
however, is the DR coupling h , which must be connected with the experimental information on 
the bottom mass: this issue will be discussed extensively in Section 4. 



3 A convenient renormalization prescription 



General low-energy analyses of the MSSM parameter space do not refer to boundary conditions 
at high scales. These analyses are usually performed in terms of parameters with a more direct 
physical interpretation, such as pole masses and appropriately defined mixing angles in the squark 
sector. Such an approach requires modifications of our two-loop formulae, Eqs. (0)^(|8|), induced 
by the variation of the one-loop parameters when moving from the DR scheme to a different 
scheme. We recall that, at the one-loop level, the two VEVs (^1,^2) and the mass parameter 
\i are not renormalized by the strong interactions. Therefore, the only parameters in the Higgs 
mass matrix that require a one-loop definition are (h , A^, S2e b ,m^, m^ ), although only four of 
these are independent, because of the relation (||). 

The sbottom masses (m~ , m~ ) in Eq. @ can be naturally identified with the pole masses. 
For the generic parameter x, we define the shift from the DR value x as 5x = x 
to this definition, we find 
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where Hij(p ) denotes the real and finite part of the (ij) component of the sbottom self-energy 
(i, j = 1, 2), and <5m~ is obtained from Eq. (||) by the interchange to- <->■ m~ . 

The most convenient definition of (/i , Ab, S2e b ) is less easily singled out. To clarify this point, 
we recall the parallel case of the 0(ata s ) corrections. In that case, besides the stop pole masses, 
the remaining independent parameters are chosen to be |4], ||] a conveniently defined stop mixing 
angle, s 2 # t , and the top Yukawa coupling h^ ole , as defined by the top pole mass M t via the 
relation Mt = h v t ole V2j ^[2. Then, the stop counterpart of Eq. (|5|) is used to establish the one- 
loop definition of At in terms of the pole top and stop masses and of the stop mixing angle. In 
the case of the 0(oeba s ) corrections, a similar procedure is not appropriate since, as can be easily 
seen from Eq. (||), S2o b is independent of Ab in the large tan/? limit. A second complication 
arises from the large one-loop threshold corrections || proportional to V2 that contribute to 
the pole bottom mass: for our calculation, the relevant ones are those 0(a s ), associated with 
one-loop SQCD diagrams with gluinos and sbottom quarks on the internal lines. As noticed in 
0, a definition of Ab in terms of the pole bottom and sbottom masses through Eq. (|J) would 
produce very large shifts in Ab with respect to its DR value, SAj, = 0(a s fj? tan 2 (5/mg). A DR 
definition for the parameters {hb, Ab, S2o b ) would avoid this problem, but would still suffer from 
the known fact that it does not make manifest the decoupling of heavy particles, for example a 
heavy gluino. 

We then look for definitions of the relevant parameters that automatically include the decou- 
pling of heavy gluinos, allow to disentangle the genuine two-loop effects from the large threshold 
corrections to the bottom mass, and provide a consistent prescription for Ab in the large tan/3 
limit. There are two quantities that have a natural physical interpretation, 

v hbV l A a \ hbV f i v hb ( a \ hbAb 

X h = -j=(cpA b + spu) — Y b = -^{spA b -Ci3fi) — ( 10 ) 

where the arrows denote as before the large tan /? limit. At the classical level, Xb is the off- 
diagonal term in the sbottom mass matrix, related to the mixing angle S2Q b via Eq. (|5|), and 
Yb is proportional to the coefficient of the trilinear {bib* R A) interaction, or, equivalently, of the 
{bib^A) interaction. 

A suitable definition of the mixing angle # , with the virtue of being infrared (IR) finite and 
gauge-independent with respect to the strong interaction, is [[Tcfl : 

I n 12 (m? )+n 12 (m? ) 
80t = -^y 5—^- , (11) 
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where IIi2(p 2 ) turns out to be independent of p 2 in the large tan/? limit. Using Eq. (^), the 
prescription on 9r can be immediately translated into a prescription for Xb- 
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Since, in the large tan (3 limit, v and are not renormalized by the strong interactions, the 
prescription on Xb can in turn be translated into a prescription for hb- Explicitly: 
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We stress that our renormalized hb, as defined above, differs at the one-loop level both from the 
DR quantity hb and from the quantity h v h ° le that would be obtained by plugging the pole bottom 
mass, M b , into the tree-level formula: 
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Concerning the definition of A b , we observe that the Yukawa coupling hb multiplying Ab can 
be absorbed in a redefinition of the trilinear soft-breaking term, A b = h^Ab. The shift in Ab 
could be defined via a physical process, e.g. one of the decays b\ — > 62 A or A — ► 61 6|, but 
such a definition would suffer from the problem of infrared (IR) singularities associated with 
gluon radiation. To overcome this problem, and given our ignorance of the MSSM spectrum, 
we find less restrictive to define 5Ab in terms of the (bib^A) proper vertex, at appropriately 
chosen external momenta and including suitable wave function corrections, so that the resulting 
combination is IR finite and gauge-independent, and gives rise to an acceptable heavy gluino 
limit. Denoting the proper vertex bibZ,A with zAi2a(Pi,P2>Pa)> we define 
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The above definition can be interpreted as the large tan (3 limit of a renormalization prescription 
on Yb, as defined in Eq. (|i~0|), since in that limit Yb — > Ab/y/2. Notice the strong resemblance with 
the corresponding renormalization prescription on X b , Eq. (|i~2|). At 0(a s ), gauge independence 
and IR finiteness follow from the fact that one-loop gluon diagrams satisfy the identity 

= iijn 11 ( P 2 )-n 2 2(p 2 )] ( 
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so that the gluon contribution to 5Ab can be written simply as 



[5 A. 



big 



A h 



U n (m 



bi- 



n 22 (m? 2 ) 



m 



61 



TO 



(17) 



b 2 



where the on-shell self-energies IIn (m| ) and I122(to? ) are indeed gauge-independent and IR 
finite. Writing 



8A b = 8h b A b + h b 5A b 



we find 
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2 This definition is suitable at (D(a s ). It can be generalized to the case of Yukawa corrections by specifying a 
prescription for the A wave function. 
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With our one-loop specifications of h b and A b , Eqs. ( |i~3| ) and (|l9|), the CP-even Higgs boson 
mass matrix takes again the form of Eqs. but the one-loop part of the corrections must 

now be evaluated in our renormalization scheme, and the functions F 21 and G 2i read now, in 
units of g 2 s C F N c / (16 vr 2 ) 2 : 
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Notice that, in this scheme, 
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in terms of our renormalized quantities (m~ ,777? , h b ,A b ), the corrections have a smooth heavy 



do not depend explicitly on Q. We also stress that, 

bi' 6 2 

gluino limit. In fact, in contrast with the case of the 0(ata s ) corrections, the gluino decouples 
for rrig — > 00, since 777 ^ G 2i — > and F 2 ^ reduces to the first line of Eq. (EDI). 



4 Input parameters 



Phenomenological analyses of the MSSM parameter space should exploit the experimental infor- 
mation on the bottom mass. Instead of expressing such information with the pole mass M b , it is 
convenient to use directly the running mass, in the SM and in the DR scheme, evaluated at the 
reference scale Qo = 175 GeV. Following a procedure outlined in [jPI]] , we take as input the SM 
bottom mass in the MS scheme, rn b (m b )^^ = 4.23 ±0.08 GeV, as determined from the T masses 
12]; we evolve it up to the scale Q$ by means of suitable renormalization group equations |Oj; 



finally, we convert it to the DR scheme. The result, which accounts for the resummation of the 
universal large QCD logarithms, is: 



m b (Q )sM 



2.74 ±0.05 GeV. 



The relation between h b = ^&(Qo)mssm anc ^ m b OI * Eq. (|22|) is given by: 
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The running parameter hb is the appropriate input quantity to be used with the DR result 
presented in Section 2, while the formulae obtained in Section 3 should be used with hb = hb~5hb, 
as defined in that section, evaluating Eq. (O) for Q = Qo- 

Notice that in Eq. (23) the large 0{a s ) threshold corrections |8j parametrized by e have 
been resummed to all orders as in [14]. With the same strategy, we can easily include the 0(at) 
threshold corrections to the bottom mass, which are expected to generate the largest two-loop 
0(atoeb) corrections to the neutral Higgs boson masses. It is sufficient to add to e the analogous 
quantity 
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where the mixing between gauginos and higgsinos has been neglected, so that the masses of the 
higgsinos coincide with /j. 

For computing the two-loop 0{ata s + at,a s ) corrected Higgs masses, as will be done in the 
numerical examples of the next section, a suitable specification must be given for the parameters 
entering the tree-level mass matrix and the one-loop 0(at + a ) corrections. In our effective 
potential approach, the tree-level mass matrix is expressed in terms of the pole mass m A and 
of the DR parameter tan/3, evaluated at the reference scale Qo, while the renormalization of 
the Z boson mass (whose numerical value we fix at m z = 91.187 GeV) does not affect the 
0(ata s + otba s ) corrections. The parameters v = 246.218 GeV and \i first appear at the one- 
loop level and do not receive corrections at 0(a s ). For the top-stop sector, we take as input the 
top pole mass, conventionally fixed at M% = 175 GeV, and the parameters (ttt-q £ ,77777, At) that 
can be derived by rotating the diagonal matrix of the On-Shell (OS) stop masses by the angle 
Of, defined as in ||. Concerning the sbottom sector, additional care is required, because of our 
non-trivial definition of hb, Eq. (|23|), and of the fact that, at 0(a s ), the parameter ttt-q ^ entering 
the sbottom mass matrix differs from the corresponding stop parameter ttiq £ by a finite shift §. 
We start by computing the renormalized coupling hb as given by Eqs. (^2|)-(|2T 



and (13). Then 

we compute nig ~ b following the prescription of || . Finally, we use the parameters hb and uIq ~ b 
to compute the actual values of the OS sbottom masses and mixing angle. The remaining input 
quantities, appearing only in the two-loop corrections, are the gluino mass rrig and the strong 
coupling constant, whose numerical value we fix at ct s (Qo) = 0.108. 
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tanP = 40 n=1.2TeV 




Figure 1: The Yukawa coupling hb, as denned in Eq. (|23|): as a function of u, for tan/3 = 40 (left 
panel); as a function of tan/3 for fi = 1.2 TeV (right panel). The other parameters are Ab = 2 
TeV, rriQ = mo = rrig = 1 TeV. The quantity h^ ole = \/2 Mb/v\ is also shown for comparison. 



5 Numerical examples 

We are now ready for some numerical examples. To prepare the ground, we study the variation 
of our renormalized hb with respect to other parameters, keeping the reference bottom mass nib 
fixed to the central value of Eq. (p^). 

The left panel of Fig. [I] shows hb as a function of fi (solid line), for tan/3 = 40. The other 
relevant parameters are chosen as Ab = 2 TeV, thq = = rrig = 1 TeV (the precise definition 
of itiq is not relevant in this case). The quantity h^ ole = ^/2Mb/v\ is also shown as a dashed 
line. The curve corresponding to hb would be very close to that of hb, thus we do not display 
it. We see that having large values of tan/3 and n, is a necessary but not sufficient condition for 
having a sizeable h^. when the threshold contribution to the bottom mass dominates, |ef,| 3> 1, 
hb must decrease for increasing values of |/i|tan/3. We also see that, when there is an almost 
complete destructive interference between the two contributions to the bottom mass, ~ — 1, 
the correct value of the bottom mass cannot be reproduced by the one-loop formula for hb in the 
perturbative regime, and the corresponding set of MSSM parameters must be discarded. Finally, 
we can see that the renormalized hb can be large only for positive ^] values of //. We then focus 
our attention on the case in which [i is large and positive, so that hb and the corresponding 
corrections to the Higgs masses can be sizeable. 

For completeness, we should mention (for recent discussions and references, see e.g. flEfi ) 
that models with b-r Yukawa coupling unification at the GUT scale favour, in our conventions, 
3 Our convention for the sign of /i is implicitly defined in Eq. (0). 
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a positive sign of Lim g , which leads to a negative e^. For sufficiently small \{j,\, radiative B decays 
and the muon anomalous magnetic moment may favour a negative sign of \x M2 , where M2 is the 
SU(2) gaugino mass, and a positive sign of fj, At. Similar but more model-dependent constraints 
can be extracted, with the help of additional assumptions on the soft supersymmetry-breaking 
terms, from the cosmological relic density. Finally, having \x and tan j3 simultaneously large may 
require a certain amount of fine-tuning [|l(| . 

The right panel of Fig. |l| shows hb as a function of tan/3, for /j, = 1.2 TeV. Again, the curve 
for hb would be practically indistinguishable and we do not show it. The other parameters are 
chosen as in the left panel, and the value of h% ole is also shown. We can see that, for this choice 
of parameters (to be taken in the following as a representative one) , values of tan (3 much larger 
than 40-50 would imply a value of hb beyond the perturbative regime. On the other hand, for 
low values of tan/3 the coupling hb is even smaller than h^ ole , and the corresponding corrections 
to the Higgs masses are expected to be negligible. For this reason, in the numerical examples of 
the 0(abCt s ) corrections we restrict ourselves to values of tan/3 between 25 and 45. 

Figs. ^ and ^ show the light Higgs mass rrih as a function of tan (3 for ix = 1.2 TeV, and as a 
function of // for tan/3 = 30, respectively. In each figure, the left panel corresponds to m A = 120 
GeV and the right panel to m A = 1 TeV . The other input parameters are chosen as At = Ab = 2 
TeV, ttiqi = rnjj = rri£> = m g = 1 TeV. For this choice of parameters, trig ^ differs from ttiqi 
by less than 1%. The curves in Figs. || and ||| correspond to the one-loop corrected rrih at 
0{a.t) (long-dashed line) and at 0(at + ab) (dot-dashed line), and to the two-loop corrected 
rrih at 0(ata s ) (short-dashed line) and at C(a t Q s + af,a s ) (solid line), respectively. We can 
see from Fig. || that, while the 0(at) prediction for rrih is practically independent of tan/3 for 
tan/3 > 25, the 0{ab) corrections lower rrih considerably when tan/3 increases. Fig. ^ shows 
that a similar decrease in rrih occurs when fi increases. Both effects are enhanced by the steep 
dependence of the renormalized coupling hb on tan /3 and /i, depicted in Fig. [I]. Comparing the 
solid and the short-dashed curves, we can see that the 'genuine' two-loop 0(aba s ) corrections 
to the Higgs mass, given by Eqs. @-@ and (|20|)-(pl|), are usually a small fraction of the O(ab) 
ones, but the former can still reach several GeV when the latter are very large. In particular, for 
small m A the 0{aba s ) corrections can be comparable in magnitude with the 0(ata s ) ones. We 
stress that the absence of very large two-loop effects from the sbottom sector is a consequence 
of our renormalization prescription, which allows to set apart the tan /3-enhanced corrections, 
resummed to all orders in the renormalized coupling hb. If we were to adopt for the sbottom 
sector the same renormalization scheme that we use for the stop sector, the dependence on tan /3 
of the one-loop corrected rrih would be smoother, but very large corrections (growing as tan 2 /?) 
would appear at two loops, questioning the validity of the perturbative expansion. 

Finally, Figs. ||a (left panel) and ||b (right panel) show both CP-even Higgs masses, rrih and 
m H , as functions of the CP-odd Higgs mass, in the region of relatively small m A (80 GeV < m A < 
180 GeV), for two different choices of the parameters. In both figures we have chosen [i = 1.2 
TeV and rriq 1 = mjj = mo = m g = 1 TeV. In Fig. ||a the other parameters are tan (3 = 40 and 
At = Ab = 2 TeV. From Fig. |^a we see that, as anticipated above, when m A is around 120 GeV 
the 0{a.bOt s ) corrections to rrih are of the same size of the 0(oeta s ) ones. This is mainly due to 

4 In the calculation of the O(at) and 0(ctb) corrections we include the effects proportional to m% and the 
momentum corrections as in ill 71] . 
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Figure 2: The mass m/, as a function of tan/3, for m A = 120 GeV (left panel) or 1 TeV (right 
panel). The other parameters are fi = 1.2 TeV, A t = = 2 TeV, wiq i = mjj = txid = m,g = 1 
TeV. The meaning of the different curves is explained in the text. 




Figure 3: The mass as a function of fx, for m A = 120 GeV (left panel) or 1 TeV (right panel). 
The other parameters are tan (3 = 30, A t = = 2 TeV, uiq i = m\j = mp = rrig = 1 TeV. The 
meaning of the different curves is explained in the text. 
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Figure 4: The masses rrih and m H as a function of m A , for /i = 1.2 TeV and mg t - = mjj = 
mo = trig = 1 TeV. The other parameters are (a) tan/3 = 40, At = Ab = 2 TeV and (b) 
tan/3 = 45, At = 1.5 TeV, ^ = 0. The meaning of the different curves is explained in the text. 

the large value of Aj,, which enhances the correction to (AA^|)5f , relevant for when m A is 
small. The 0(ab<y s ) corrections to m H are rather small in this example, but they can be larger 
for different parameter choices. In Fig. ||b the relevant parameters are tan/3 = 45 , At = 1.5 
TeV and Ab = 0. For this choice, radiative corrections mainly affect (AA1|)^. Thus one of 
the eigenvalues is roughly degenerate with m A and receives small corrections, while the other 
eigenvalue is almost independent of m A and receives large corrections. In particular, the genuine 
0(abCt s ) corrections to either or m H are around 3 GeV in this example. 

6 Conclusions and discussion 

In this paper we presented explicit and general results for the 0{a.b a s ) corrections to the MSSM 
neutral Higgs boson masses, in the physically relevant limit of large tan (3. Actually, a large 
value of tan (3 is a necessary but not a sufficient condition for having large corrections, which 
require sizeable values of both \x and hb- We proposed a renormalization prescription for the 
sbottom sector that automatically includes the decoupling of heavy gluinos and separates the 
large threshold corrections, appearing in the relation between hb and the pole bottom mass, from 
the genuine two-loop effects. We also discussed the numerical impact of our results in a number 
of representative examples. 

A complete study of the two-loop (s)bottom corrections would require also the knowledge of 
the 0(a.bOit) and O(o^) effects. Concerning the former, it is plausible that the most important 
effects can be taken into account by adding to e& the analogous quantity e fc '. The 0(cef ) corrections 
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would need a dedicated calculation, but an estimate of their importance can be obtained from our 
knowledge of the 0(a 2 ) corrections. In Refs. ||, [7||, explicit formulae for the 0(a 2 ) corrections to 
the Higgs masses, valid under simplifying assumptions on the MSSM parameters, were presented. 
The corresponding formulae for the 0(a 2 ) corrections can be derived from such formulae by 
performing suitable substitutions and taking appropriate limits. In the case of large tan/3 and 
universal soft sbottom masses, degenerate with m A and much larger than the weak scale (mg = 
mo = m A = Ms S> v), it is possible to derive a simple expression for the Oia^ + a s + a 2 ) 
corrections to m?: 



4a. s 



2 ' 

mi 



8tt 2 M s v 2 



1 + ^fnjA) + 



3 a.b 

47T 



f \M 2 



A 2 



2 ln^J- + - ln2 + 1 | +C , 

v Ml 3 M 2 



(27) 



where = h^v fi/V2 in the large tan (3 limit, C ~ 0.27 and f(x) is a positive function, defined 



as 



x(3x-2) x 2 (3x 2 -8x + 6) 



x — 1 



1) 



Inx + (3 x + 1) (x — 1) In |x — 1| 



(28) 



Some limiting values are f(0) = 0, f(l) = 9/2, /(oo) = 3/2. In view of the result in Eq. (J27 
we expect that, for values of a& not much larger than a s , the 0{a 2 ) corrections should be at 
most comparable with the 'genuine' 0(ai,a s ) effects. 



Acknowledgments 

P.S. thanks A. Dedes, A. Quadt and V. Spanos for discussions. F.Z. thanks E. Franco and 
G. Martinelli for discussions, the Physics Department of the University of Padua for its hospitality 
during part of this project, and INFN, Sezione di Padova, for partial travel support. This work 
was partially supported by the European Programmes HPRN-CT-2000-00149 (Collider Physics) 
and HPRN-CT-2000-00148 (Across the Energy Frontier). 



References 

[1] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, The Higgs Hunter's Guide, Addison- 
Wesley, 1990 and (errata) |hep-ph/9302"272 . 

[2] J. Ellis, G. Ridolfi and F. Zwirner, Phys. Lett. B257 (1991) 83 and Phys. Lett. B262 (1991) 
477; 

Y. Okada, M. Yamaguchi and T. Yanagida, Prog. Theor. Phys. 85 (1991) 1 and Phys. Lett. 
B262 (1991) 54; 

H. E. Haber and R. Hempfling, Phys. Rev. Lett. 66 (1991) 1815. 
[3] R. Hempfling and A. H. Hoang, Phys. Lett. B331 (1994) 99 ||hep-ph/94012Tg ] . 



13 



[4] S. Heinemeyer, W. Hollik and G. Weiglein, Phys. Rev. D58 (1998) 091701 [|hep-ph/9803277[| 
Phys. Lett. B440 (1998) 296 [|hep-ph/9807423f , Eur. Phys. J. C9 (1999) 343 jhip~- 
ph/9812472] , and Phys. Lett. B455 (1999) 179 [ |hep-ph/ 9903404] ]; 



R. Zhang, Phys. Lett. B447 (1999) 89 |frep-ph/9808299| ; 

J. R. Espinosa and R. Zhang, JHEP 0003 (2000) 026 jhep-ph/991223l| . 



[5] G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B611 (2001) 403 |hep-ph/0105096(| . 
[6] J. R. Espinosa and R. Zhang, Nucl. Phys. B586 (2000) 3 jhep-ph/0003246[| . 



[7] A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B631 (2002) 195 jhep- 



ph/0112177|1 . 



[8] T. Banks, Nucl. Phys. B303 (1988) 172; 

L. J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D50 (1994) 7048 [|hep-ph/9306309 
R. Hempfling, Phys. Rev. D49 (1994) 6168; 

M. Carena, M. Olechowski, S. Pokorski and C. E. Wagner, Nucl. Phys. B426 (1994) 269 
l |hep-ph/9402253 |. 

[9] A. Bartl, H. Eberl, K. Hidaka, T. Kon, W. Majerotto and Y. Yamada, Phys. Lett. B402 
(1997) 303 [|hep-ph/9701398l ; 

H. Eberl, K. Hidaka, S. Kraml, W. Majerotto and Y. Yamada, Phys. Rev. D62 (2000) 055006 
|iep-ph/9912463 |. 

[10] A. Pilaftsis, Nucl. Phys. B504 (1997) 61 jhep-ph/9702393(1 ; 



J. Guasch, J. Sola and W. Hollik, Phys. Lett. B437 (1998) 88 [|hep-ph/9802329|| ; 
H. Eberl, S. Kraml and W. Majerotto, JHEP 9905 (1999) 016 |hep-ph/990341~l ; 
Y. Yamada, Phys. Rev. D64 (2001) 036008 [|hep-ph/010304^ . 

[11] J. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Srednicki, Phys. Lett. B510 (2001) 236 
[|hep-ph/0102098 |. 



[12] M. Beneke and A. Signer, Phys. Lett. B471 (1999) 233 flhep-ph/9906475| ; 

A. Hoang, Phys. Rev. D61 (2000) 034005 jhep-ph/990555CH and |iep-ph/0008l"02 



[13] K. G. Chetyrkin, Phys. Lett. B404 (1997) 161 [|hep-ph/9703278 |. 



[14] M. Carena, D. Garcia, U. Nierste and C. E. Wagner, Nucl. Phys. B577 (2000) 88 jrep- 



| ph/99125iq] ; 

G. Degrassi, P. Gambino and G. F. Giudice, JHEP 0012 (2000) 009 ||hep-ph/0009337| . 

[15] U. Chattopadhyay, A. Corsetti and P. Nath, |hep-ph/02042~5l"| ; 

H. Baer, C. Balazs, A. Belyaev, J. K. Mizukoshi, X. Tata and Y. Wang, JHEP 0207 (2002) 
050 Hhep-ph/0205325 1. 



[16] A. E. Nelson and L. Randall, Phys. Lett. B316 (1993) 516 jhep-ph/9308277 [. 
[17] A. Brignole, Phys. Lett. B281 (1992) 284. 



14 



